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The zero momentum sectors in effective theories of three dimensional QCD coupled to pseudo- 
reaf (two colors) and real (adjoint) quarks in a classically parity-invariant manner have alternative 
descriptions in terms of orthogonal and symplectic ensembles of random matrices. Using this cor- 
respondence, we compute finite- volume QCD partition functions and correlation functions of Dirac 
operator eigenvalues in a presence of finite quark masses of the order of the smallest Dirac eigenvalue. 
These novel correlation functions, expressed in terms of quaternion determinants, are reduced to 
conventional results for the Gaussian ensembles in the quenched limit. 
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I. INTRODUCTION 

Spontaneous breaking of global symmetries has long 
been a subject of an extensive study for quantum field 
theories in various dimensions. Besides well-known cases 
such as 2D Gross-Neveu and Schwinger models and 4D 
QCD that exhibit breakdown of discrete or continuous 
chiral symmetry, it has been suspected that 3D QED 
and QCD may also undergo spontaneous breaking of fla- 
vor symmetry (l]-Q. For such odd dimensional field the- 
ories, there is no comprehensive theorem that predicts 
the surviving part of the flavor symmetry as in the case 
of even dimensions || , as the notion of fermion chirality 
is absent. However, if there are even Nf — 2n num- 
ber of massless 2-component complex fermions ipj, one 
can group them appropriately into 4-componcnt complex 
fermions * 3 [§-[§: 



C = 



In 

i=i 



2*i^7p*i, (1-1) 



where 7 M = ay, ® 03 (/i = 1,2,3) in Euclidean 3D 
space. Then one can define a gwasi-chirality accord- 
ing to two Hermitian matrices, e.g. 74 = 1 (S> ci and 
75 = 1 g) (72, that anticommute with the Dirac opera- 
tor -0^7^ and with each other. Generators of the flavor 
U(2n) group are rearranged into that ofU(n) group times 
{1,74,75,7s = 1 ® 03} in the 4-component notation. In 
order to predict spontaneous flavor symmetry breaking 
pattern along the same line as in 4D QCD, one introduces 
a small symmetry-breaking mass term that is parity- 
invariant, i.e., ^2j—i mj^fj^j but not Xw=i mj^jjs^j- 
In the 2-component notation it leads to include masses 
{to} = (mi, . . . , to„, —mi, . . . , —ran). As the fermion de- 
terminant 



J|dct(iL>, 



7m 



= Y[c\ct(-(D n ,cr f ,) 2 +mj) (1.2) 



is positive definite under this restriction, one can appeal 
to the Vafa-Witten theorem fjl and predict that if the 
flavor symmetry is spontaneously broken, the absolute 
values of the fermion condensate (ipj^pj) are equal for all 
j = 1, . . . , 2n and their signs are the same as those of 
respective masses. That is, the continuous part of the 
global symmetry group is broken according to [Hp) 



U (2n) -> U(n) x U{n) 



by the order parameter 



1 In 



(1.3) 



(1.4) 



while the discrete Z2 group (a product of the parity and 
the exchange of fields tpj <-> ip n +j) remains unbroken. 
The formation of the quark condensate was indeed ob- 
served in Monte-Carlo simulations on a lattice 0. This 
peculiar pattern of flavor symmetry breaking can also be 
predicted for 3D large- N c QCD by the Coleman- Witten 
argument B , as remarked in Ref. |6j . One can repeat the 
above argument valid for complex (fundamental repre- 
sentation of SU (N c > 3) gauge group) fermions towards 
the cases with even flavors of pseudoreal (fundamental 
representation of SU(2) gauge group) and of real (ad- 
joint representation of SU(N C ) gauge group) fermions. 
We assign Dyson indices [3 = 2, 1, 4, respectively, to these 
three cases, according to the anti-unitary symmetries of 
the associated Dirac operators M. Then the continuous 
parts of the global symmetry groups are predicted to be 
broken down as [11011 



Sp(2n) 
SO(2n) 



Sp(n) x Sp(n) 
SO{n) x SO(n) 



09 = 1), 
09 = 4). 



(1.5a) 
(1.5b) 



These symmetry breaking patterns determine the forms 
of the low-energy effective Lagrangian of associated 
Nambu-Goldstonc bosons: 



1 



C cS = — tr dfJJd^ - Etrmt/rC/ 1 + ■••, (1.6) 

where U(x) takes its value in the coset manifolds^ 

M F = AIII„, n , CII n ,„, BDI n ,„, (1.7) 

respectively for f3 = 2, 1, 4. The mass and chiral matrices 
in the above are defined as 



m = 



diag(mi, . 
diag(mi, . 

In 
In 



03 
0- 3 



. ,m n ) 
. ,m n ) 

ifi 
(fi 



2,4) 
1) 



ia 2 



03 
09 



2,4) 
1) 



(1.8) 
(1.9) 



Given non-linear a models (NLcrMs) in 4D, Leutwyler 
and Smilga fll|| proposed to extract out of them nonper- 
turbative, exact information on the spectra of Dirac oper- 
ators by imposing a constraint on the parameters: the lin- 
ear dimension L of the system be much smaller than the 
Compton length of Nambu-Goldstone bosons ~ F/ VmE. 
In this 'ergodic' regime where the zero momentum mode 
of U (x) dominates, the effective finite- volume partition 
function simplifies into a finite-dimensional integral and 
its small-m expansion yields a sequence of spectral sum 
rules. Verbaarschot and collaborators [|l2|]l5| l (see Ref. 
fuj for a exhaustive list of references) made an impor- 
tant observation along this line, that the OD finite- volume 
partition functions could as well be derived from mod- 
els much simpler than QCD, random matrix ensembles 
(RMEs). In our 3D QCD context it means that 

in the limit 

L — > do, rrii — > 0, fa = L 3 Em; : fixed, (1-10) 
the finite-volume partition functions 

Z{{^})= ( dU exp(tr/xC/rt/ t ) , (1.11) 

JMf 

{i = L 3 Em, 



have alternative representation in terms of large- N RMEs 

« n 

Z({m})= dHe- puv ^T\&et(H 2 + mJ), (1.12) 

where the integral domains V are sets of N x N complex 
hermitian (u(N) = T(Ajv)), real symmetric (o(N) — 
T(AIjv)), and quaternion selfdual (sp(N) = T(AIIjv)) 
matrices H for (3 = 2, 1, 4, respectively. The determinant 
in the (j = 4 case is understood as a quaternion determi- 
nant (Tdet). Their proofs consist of the 'color- flavor' (or 
Hubbard- Stratonovich) transformation p5|Jl^ ] that con- 
verts the integration variables into matrices with small 
(n x n) dimensions, and the saddle point method under 
which 

N — * oo, rrii — > 0, fa = irp(0)mi : fixed. (1.13) 



Here p(0) stands for the large- TV spectral density of the 
random matrix H : 



p{x) = lim 

N —>o> 



(trS(x-H)), 



(1.14) 



at the spectral origin. The RMEs (1.12) are motivated by 
the microscopic theories (3D Euclidean QCD) on a lat- 
tice, with a crude simplification of replacing matrix ele- 
ments of the Hermitian Dirac operator ip = (id^+A^a^ 
by random numbers Hjk distributed according to the 
weight e~P trV ( H ). U nder this correspondence, the mi- 



crosco pic l imit (1.13) is equivalent to Leutwyler-Smilga 
limit fll.lC| ), since the size N of the matrix H is inter- 
preted as the number of sites L 3 of the lattice on which 
QCD is discretized, and the Dirac spectral density at zero 
virtuality p(0) is related to the quark condensate by the 
Banks-Casher relation S = irp(0)/L 3 @. 

On the other hand, the situation is far more subtle 
in the case with an odd number (Nt = In + 1) of 2- 
component fermion flavors l^fl. In the massless case, 
the fermion determinant det ™ +1 (iD Al a„) is not posi- 
tive definite, and its phase is gauge dependent. This 
dependence can be compensated by the Chern-Simons 
term that is yielded by a gauge-invariant regulariza- 
tion, but this anomalous term explicitly breaks the par- 
ity B. Therefore, even in a presence of small masses 
{m} = (mi, . . . , m n , —mi, . . . , — m„, 0) that respects the 
2/2 invariance (combining parity and flavor exchange) 
classically, one cannot appeal to the previous argument 
to derive low-energy effective Lagrangians. With this 
situation in mind, we nevertheless adopt a pattern of 
spontaneous flavor symmetry breaking proposed by Ver- 
baarschot and Zahed || for (3 = 2 and its generalizations 
to /3 = 1 and 4: 

U(2n + 1) -»■ U{n) X U{n + 1) (0 = 2), (1.15a) 
Sp(2n + 1) -> Sp(n) x Sp{n + 1) 03=1), (1.15b) 
SO(2n + l) -» SO(n) x SO(n + 1) [fi = A), (1.15c) 

leading to NLcrMs of Nambu-Goldstone fields over the 
coset manifolds 



M, 



AIII 



n,n+l i 



CII 



i, BDI nn _(_i, 



(1.16) 



respectively. Then, by the same token as in the case of 
even Nf, one can write down corresponding RMEs [fl8| : 

„ n 
Z({m})= / dHe^ trV ^ det TT det (if 2 + m|) , 



i=l 



(1.17) 



which are equivalent, in the limit (1.13), to the 'finite- 
volume partition functions' Iq] 



2(M) = 



dU 



Mi 



cosh(tr/x£/r[/t) (JV : even 
sinh(tr^UTU't) (N : odd) 



, (1.18) 



. diag(l„,-l„+i) {fi = 2,4) , , 

< diag(l„,-l n+ i)«8>i<7 2 08 = 1) • [L - W> 
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As the Chern-Simons term cannot be incorporated within 
these RMEs, their physical relevance is unclear [[|. An 
immediate problem is that if the ran k N of th e ma - 
trix H is odd, the partition function ( 1.17 ) (or ( |l.l<j )) 
is zero, and the (unnormalized) correlation functions 
are odd under a simultaneous change of signs of the 
arguments, which are unacceptable as physical observ- 
ables. The above relationships between RMEs and 
NLoMs for 3D QCD consist, together with its coun- 
terpart for 4D QCD (M F = A„,AII„,AI„ and V = 
T{Aill N , N ,),T{Vm N . N ,),T(Cll N , N ,) for (3 = 2, 1,4, re- 
spectively), a part of Zirnbauer's complete classification 
scheme of RMEs in terms of Riemannian symmetric 
spaces pJ| ]. 

These RMEs are technically suited for the computation 
of correlations of eigenvalues {x} of the Dirac operator 
ip r*> H in the microscopic asymptotic limit where the 
ener gy ei genvalues are scaled as the quark masses are in 

Eq.fri, 



N 



co, x — * 0, A = irp(0)x : fixed. 



(1.20) 



For the chiral RMEs describing 4D QCD in the low- 
energy ergodic regime, such Dirac spectral correla- 
tors have been computed previously for the massless 
jl3],^0-24| and recently for the massive cases p5|-|3l|| with 
all three values of (3. On the other hand, for the (non- 
chiral) RMEs describing 3D QCD, they have been analyt- 
ically computed solely for the unitary ((3 = 2) ensemble, 
in the massless case [p0p|] , as well as in the massive case 
p2| , fL8f . For other values of (3, a numerical work based on 
a finite- N formula for the correlation functions recently 
appeared only in the massless case Q. The subject of 
this Article is to complete this program by analytically 
computing the partition and correlation functions for the 



orthogonal ((3 = 1) and symplectic (f3 = 4) ensembles in 
a presence of finite scaled mass parameters. We employ 
a slightly modified version of the method used in our 
previous articles [^9| , ^0| . 

We finally remark on the universality issue. It was 
noticed by §ener and Verbaarschot |34j (see also Ref. 
|p5| ) and proved by Widom |3(J that the diagonal ele- 
ment S(x, y) of the quaternion kernel for an orthogonal 
or symplectic ensemble, and accordingly all spectral cor- 
relation functions thereof, can be constructed from the 
scalar kernel K(x, y) for a unitary ensemble sharing the 
same weight function: 



S T = (i — (I — K)sKDJ K (/3=1), (1.21a) 
S = (i - (I - K)DKe) K (/? = 4), (1.21b) 



where /, D, e, S, K stand for integral operators with 
convolution kernels S(x — y), 8'(x — y), ^sgn(x — y), 
S(x,y), K(x,y), respectively, T and _1 stand for trans- 
pose and inverse opera to rs. S ince the scalar kernel in the 
asymptotic limit ( |l.20| ), ( 1.13 ) is insensitive to the details 
of the potential V(x 2 ) either in the absence |57],[58| or in 
the presence of finite and nonzero /i's [ ^2|Jl8| , the univer- 
sality of correlation functions for orthogonal and sym- 
plectic ensembles are automatically guaranteed"' 1 . There- 
fore it suffices for us to concentrate onto Gaussian en- 
sembles, V(x 2 ) = x 2 /2. This choice leads to Wigner's 
semi-circle law 



p(x) = - \/2N -x 2 . 



(1.22) 



II. ORTHOGONAL ENSEMBLE 



For (3 = 1, we treat the following three cases separately: 



A 
B 

C 



{m} = (mi, . . . , m„, -mi, 
{m} = (mi, . . . , m„, -mi, 
{m} = (mi, . . . ,m„, -mi, 



-m n ), 

-m„, 0), even N, 
-m n , 0), odd N. 



A. even Nf 



We first consider the case with N f = 2 n flavors and {m} = (mi, . . . , m„, —mi, . . . , — m n ). 

We express the partition function ( 1.1 2| ) of the RME in terms of eigenvalues {xj} of H (up to a constant independent 
of m): 



-i roo poo N N , n \ N 

■kf-f n^nH /2 n^+^))ni^ 

J -°°J = 1 3 = 1 V fc=l 7 j>k 



Z({m}) 

The p-level correlation function of the matrix H is defined as 



Xk\- 



(2.1) 



3 



p(xi,...,x p ;{m}) = (~Y[tr 8(xj - H)) 

3 = 1 

_ E p (xi,...,x p ;{m}) 



E ({m}) 



E p (xi,...,x p ;{m}) = 



1 



(N-p)\ 

(So = Z). We define new variables 





[' 00 


N 


N 








n 


dxj Y[ 


(•■ 


J — oo 


3 


=p+i 


3=1 




Z2j-\ 


= irrij 




= 1,... 


n), 


Z 2 j = 


—irrij 




= 1,... 


,n), 


Z 2 n+j 


= Xj 


u 


= 1,... 


,P). 



N 



tj Xk I 



(2.2) 
(2.3) 



(2.4) 



Then the multiple integral ( |2.3| ) is expressed as 

1 



Ujll \fw{zj) Ujlk( Z 3 - Z k) 
roc poo 2n+W 



(iV- 



^ ,.00 ,.00 2n+N 2n+JV 2n+Af 2n+7V 

— — / ••• / n dz 3 n y™^) n o^-^) n ^o^-^)' ( 2 - 5 ) 

Pl- J-00 J-oOj =2n+p+1 J = 1 - >k j>k>2n 



where w(z) — e~ z . Eq.(2.5) resembles an (2n + p)-level correlation function of the conventional Gaussian ensemble 
with 2n + N levels. However, conventionally the levels zt, . . . , Z2n+p are all real, while in the present case some of 
them (z\, ... , Z2n) are pure imaginary. We carefully incorporate this fact into the following evaluation. 
Let us denote the integrand in Eq.(|2.5[) as 



2n+N 



_2n+N 



2n+N 



p(z u z 2n+N ) = n v w(yz ^ n & ~ ^ n sgn ( z -? ~ 

j—l j>k j>k>2n 



For N even, an identity 



2n+N 



Yl sga(zj - Zk) = Pf[sgn(z fe - Zj)]j,k=i,...,2n+N 
j>k 

holds for real Zi, ... , Z2 n +N- By taking the limit z\ < Z2 < ■ ■ ■ < Z2 n — * —00, we find another identity 

2n+N 

Yl Sgn(z i - Z k ) = ~Pf[Fjk]j,k=l,...,2n+N, 



(2.6) 



(2.7) 



(2.8) 



j>k>2n 



where 



sgn(fc-j) (j,k = l,...,2n) 



F 



jk 



1 



(j = 1, . . . , 2n\ k = 2n + 1, . . . , 2n + N) 
(j = 2n + 1, . . . , 2n + N; k = 1, . . . , 2n) 



(2.9) 



sgn(z fe - zj) (j, k = 2n + 1, . . . , 2n + N) 



Substitution of Eq.fl2.8D in to Eq.(|2.6D yields 



2n+N 



_2n+N 



p(zi, . . . , Z 2u +n) = Yl V W ( Z ^ II ( Z 3 _ Z fc) Pf [Fjk]j,k=l,...,2n+N- 

3=1 j>k 



For N odd, we similarly obtain 



2n+N 



P(ZX, . . . ,Z 2n+N ) = Y[ \J W ( Z 3) II ( Z 3 - Z fe) Pf 
j=l j>k 



_2n+N 



[Fjk]j,k=l,...,2n+N [gj]j=l,...,2n+N 
[-fl , fe] fc =i,...,2n+JV 



(2.10) 



(2-11) 
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with Q j = Qk = 1 (J, k = 1, . . . , 2n + N). The Pfaffians in the above can be represented as quaternion determinants 
||,||]39[Q. In doing so, we introduce monic skew-orthogonal polynomials Rj(z) — z 3 + ■ ■ ■, which satisfy the 
skew-orthogonality relation: 



(i?2j, R2k+1)R = — {Rik+l: R2j)n = r j5jk, 
(R2j,R2k/R = (Ifaj+1, R2k+1/R = 0, 



where 



{f,9)R= I dzy / w(z)g(z) / dz'y/w(z')f(z') -(/<->• g). 



(2.12) 



(2.13) 



Explicit forms for the skew-orthogonal polynomials and their norms associated with the Gaussian weight w(z) are 
known Efl]: 



R2j[z) = -2jH 2 j{z), 



#2 i+ i(z) = ^(H 2j+1 {z) - H' 2j {z)), 



= 2-^ +1 (2j)!v^, 



in terms of the Hermite polynomials 

Now we present the following theorems: 
Theorem 1 

For even N, we can rewrite p(z%, . . . , Z2ti+n) as 



= (-l)V 



dz 3 ' 



n+N/2-1 

p(z±, . . . , Z 2n +N) = ( Yl r j^ T ^ et [fjk( Z j,Z k )]j.k=l,...,2n+N- 
3=0 



The quaternion elements fjk{zj,Zk) are represented as 

fjk(zj,z k ) 



S(zj,Zk) I{zj,Zk) 
D(zj,z k ) S(z k ,Zj) 



The functions S(zj, Zk), D(zj, Zk) and I(zj, Zu) are given by 

n+JV/2-1 



S(z 3 ,z k ) 
D(zj,Zk) 

I(Zj,Z k ) 



E 

n+TV/2-1 



^(^)*2^+i(^fe) - ^2e+i{zj)^2e(zk) 



rt 



E 



*2€(«i)*M+l(«fc) - ^2t+l{Zj)^2l{Zk) 



1=0 

n+N /2-1 



'7 



E 

£=0 



^feO^+l^fc) - $2«+l(Zj)$2^(.Zfc) 



'V: 



3*' 



where 



9j(z) = ^Mz)R 3 {z), 

dz^Jw(z)R (z) sgn(z fe — z) (k = 2n + 1, . . . , 2n + N) 
— J_ dzy / w(z)Rj(z) (k : otherwise) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 
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Theorem 2 

For odd N, we have 



n+[N/2]-l 



p(Zl,...,Z 2n+N ) = ( Y[ r j ^S2n+N-lTdet[f°k d {z jl Z k )}j,k=l,...,2n+N- 



3=0 



The quaternion elements are represented as 



fT{z^z k ) = 



and 



S oM {z h z k ) I° dd (z 3 ,z k ) 
D° dd ( Zj ,z k ) S odd {z k ,z 3 ) 



^j(z)dz. 



The functions S° , D and I are given in terms of 5, D and / defined in Eq.( 2.18 ), according to 

^2n+N-l(Zk) 



S odd (z J1 z k ) = S(z J ,z k ) 



D odd {z 3 ,z k ) = D{z^z k ) 



S2n+N-1 



r aa { Zjl z k ) = i{ Zj ,z k ) 

Here * stands for a set of substitutions 



$2n+N-l{Zj) - §2n+N-l(z k ) 
S2n+N-1 



R j (z)^R j (z) 

S2n+N-1 

associated with a change in the upper limit of the sum 

N 1 
n + — l i ► n ■ 



-R 2n+N -x(z) (j=0,...,2n + JV-2), 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



Theorem 3 

Let the quaternion elements qj k of a selfdual n x n matrix Q„ depend on n real or complex variables z±, 

Qjk = f jk (z j7 z k ). 
We assume that fj k (zj,z k ) satisfies the following conditions. 

fnn(z n , Z n ^dfl(z n ^ C n , 

f jn (zj,z n )f nk (z n , z k )d^{z n ) = fjk{zj,Zk) + Xfjk(zj,Zk) - fj k (zj,z k )X. 
Here dfi(z) is a suitable measure, c„ is a constant scalar, and A is a constant quaternion. Then we have 

J Tdet Q n dfi(z n ) = (c n — n + l)TdetQ n _i, 
where Q n -i is the (n — 1) x (n — 1) matrix obtained by removing the row and the column which contain z n . 



, Zn as 
(2.26) 

(2.27a) 
(2.27b) 

(2.28) 



It is straightforward to show that the quaternion elements fj k {zj, z k) and f° k {zj, z k ) in Theorem 1 and Theorem 
2, respectively, both satisfy the conditions imposed on fj k (zj,z k ) in Theorem 3. This means that we can inductively 
write 



p+[JV/2]-l 



Sp(zi, . , Z2n+p) 



11.7=0 



Tdet[/ ife (z i ,z fc )Ufe=i,...,2n+p (iV : even) 

odd^. „. M . , . _ . (AT . „AA\ ■ 



n 2 =i Vw{^)U 2 lk( z J ~ z k) I S2 «+^-i Tdet [/j 9 fc (^.■Zfc)]j,fc=l,...,2n+p (JV : Odd) 



Since the final result in the asymptotic limit N — > oo sho uld be insensitive to the parity of -/V, we consider only even 
N henceforth. Then the p-level correlation function (2.2) is written as 



,x p ;{m}) 



Ep(zi,. . . , z 2n + P ) _ Tdet[f jk (zj,Zk)]j,k= 



•l,...,2n+p 



E (zi, . . .,z 2n ) Tdet[fjk(zj,Zk)]j,k=i,...,2n 



(2.30) 



We introduce a set of notations 



Sjk — S(zj,Zk) 
Sjk — S(zj, z 2n +k) 



Sjk — S{z 2n +j, z k ) 
Sjl R — & (z 2n +j , z 2 



n+k, 



(j,k = l,...,2n), 
(j = 1, . . .,2n; k = 1, ... ,p), 
(j = l,...,p; k=l,.. . ,2tt,), 
(j,k = l,...,p) 



(2.31) 



and similarly for D and /. Using Dyson's identity (JA. 6|) we can rewrite the correlation functions as 



Pf 



p(x 1 ,...,x p ;{m}) = (-ir^ 2 - 



-I 11 S 11 —i IR s IR 

-(S H ) T D 11 -{S RI ) T D IR 

—I RI S RI —j RR s RR 

-(S IR ) T D RI ~{S RR ) T D RR 



Pf 



Pf 



= (_l)P(p-l)/2_ 



-I 11 s u 
-(S H ) T D H 

D" -{S RI ) T -(D RI ) T 



S RJ 
D RI 



_jRR S RR 
-(S RR ) T D RR 



Pf [D 11 ] 



(2.32) 



In the last line we have exploited a Pfafnan identity that holds for antisymmetric matrices A, B of even ranks and a 
row vector v: 



Pf 





V 




A 








V 




—v 1 ■ ■ ■ — v 1 


B 





PfL4] Pf[B]. 



(2.33) 



Now we proceed to evaluate the component functions of the quaternion kernel fjk{zj,z k ) in the asymptotic limit 
(1.13) and (1.20), where unfolded microscopic variables 



<2Nz 



2j-l 



>2Nz- 



2j = (,2j 



<2Nz 2n+ j — A 



(.1 = 


1,. 


. ,n), 


u = 


1,. 


.,n), 


(j = 


1,- 


■,P) 



(2.34) 



are kept fixed. We note that all elements of the sub-matrices that appear in the second line of Eq.(2.32) are expressed 
as (derivatives or integrals of) an analytic function 



n+JV/2-l 



S(z,z') = J2 

3=0 



-(, 2 + z ' 2 )/2 H 2n+N (z)H 2n+N ^(z') - H 2n+N ^(z)H 2n+N (z') 



2 2n+N ^T(2n + N) 

e^'V 2 
+ 2 2n + N ^T(2n + N)' 



H 2n+ N-i(z) / e u /2 H 2n+N (u)du, 



(2.35) 
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where 



tyj(u)du, 



(2.36) 



with real and/or imaginary arguments: 



•,11 



qRI _ 
S jk - 

D% = 

qRR _ 
&jk - 

jRR _ 
1 jk ~ 



1 9 S( \ 

7C ^ \Zj i Z k J ; 

2 w 


(2.37a) 


S(xj , Zfc), 


(2.37b) 


2 ax/ ^' Zfej ' 


(2.37c) 




(2.37d) 


1 9 ^ 

-— b[Xj ,X k ), 


(2.37e) 


—2 / S(xj, x)dx — sgn(xj — x k ). 

J Xj 


(2.37f) 



In the second line of Eq. (2.35), we have singled out the unitary scalar kernel and applied to it the Christoffel-Darboux 
formula. Substituting asymptotic formulas for the Hermite polynomials 



H 2k (z) ~ 1 ; _ ' cos(2v^z), 



H 2k +i{z) 



V irk 
{-l) k 2 2k+1 k\ . 



.(2Vkz), 



(2.38) 



valid under k — + oo, z — + 0, ykz : fixed, one can show that S(z, z') approaches the sine kernel 



/2N W2N V2N 7r(C — CO 



(2.39) 



The Pfaffian elements in the asymptotic limit, after taking into account an unfolding by the factor \/2N , are then 
expressed in terms of K((): 



D 



11 



5 J* 



D 



RI 

jk 



cRR 



D 



RR 

jk 



lf k R = I( Xj ,x k ) ~ 2 / K{\)d\ - sgn(A, - A fe ). 



^D(z 3 ,z k ) ~ ^'«? - Cfe), 



1 



/27V 



S(xj,Zk) ~ #(Aj - Cfe), 



-^D( Xj ,z k ) ~ 2^'(Aj - Cfc), 



1 



/27V 



S(xj,x k ) ~ K(Xj - A fe ), 



^( W fe)~ ^'(A, -Afe) 



(2.40a) 
(2.40b) 

(2.40c) 
(2.40d) 

(2.40e) 

(2.40f) 



These matrix elements constitute the finite- volume partition function 



Z(M)eeS ({-^=}) 



const. 



Pf[D 



7/1 



(2.41) 



and the scaled spectral correlation functions 
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/ 1 \p 
p s (Ai, . . . , A p ; {/i}) = ) P( 



Pf 



p(p-l)/2_ 





A p n 


}) 




v/i/V' L \/2iV 




-(S fl/ ) T - 


(D RI ) T ' 




_jRR 


cRR 




-(S RR ) T 


D RR 


Pf [D"l 



(-1) 



oo when the ratio of two Pfaffians is replaced by a minor Pf 



(2.42) 



In the quenched limit hi, . . . ,/x„ 

the correlation functions approach those of the Gaussian orthogonal ensemble 
sequence 



jrtrt g JTL ft 

-(S RR ) T D« H 



•,-J"ti) 



p s ({A};/xi, 



, Mn-l, — Ml, 



By the same token, it satisfies a 

(2.43) 



'Pn-l 



as each of the masses is decoupled by going to infinity. To illustrate this decoupling, we exhibit in FIG.l a plot of the 
spectral density p s (A; /i, — /i) (p = 1, n = 1). 



B. odd N f , even N 



Next we consider the case with Nf = In + 1 flavors, {to} = (toi, . . . , m n , — TOi, . . . , —m n , 0), and with even JV. We 
express the partition function (1.17) of the RME in terms of eigenvalues {xj} of H: 



Z({m}) 



N 



oo p oo JV 

■••/ n^n 

-oo ^-oo 3 = 1 j=1 



e x a/ 2 a;. 



fc=l 



j>k 



The p-level correlation function of the matrix H is defined as 

p 

p(x 1 ,...,x p ;{m}) = (JJtrtffo -if)) 



a p (xi, . . . ,x p ; {m}) 
S ({m}) 



Sp(a;i,...,a:p;{m}) 



(S = Z). We define new variables zj as 



bv/-7 n ^nk^n( 

P)-J-oo J-oo j=p+1 j=1 \ k=1 



2 . 2- 
X ■ + TO fe , 



JV 

m 

j>k 



Xj %k\ 



Zq = 0, 










U 


= 1,.. 


. ,n) 


Z2j = —irrij 


u 


= 1,. 


.,n) 


%2n+j Xj 


u 


= 1,. 


;P) 



Then the multiple integral ( 2.46| ) is expressed as 

1 



^p( z 0, ■ ■ ■ , ^2n+p) 



IL=o y/ w ( z i) Yljlk>o( z i - z k) 



r x poo 2n+N 



2n+N 

n dz i 

■•- J °° j=2n+p+l j=0 



2n+N 



2n+N 



| pOO poo ^.t-r^v ^tt-rn ^'*TJ. i«TJ. 

fv^v/ "/ n v w fe) n n sg^-**) 

V JV jCy- J —OO J —OO __n„ i „ i 1 __n i.^r. 7.^ o„ 



j>k>0 



j>k>2n 



where w(z) = c z . 

Let us denote the integrand in Eq.( [2.48 ) as 



2n+JV 



. 2n+JV 



2n+7V 



p{Z0,...,Zl n+N ) = Yl \ W ( Z 3) II ( Z 3 - Z k) Yl S S n ( Z j- Z k)- 
j=0 j>k>0 j>k>2n 



(2.44) 



(2.45) 
(2.46) 



(2.47) 



(2.48) 



(2.49) 
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An identity 



2n+N+l 

Yl sgn(zj - z k ) = Pf[sgn(z fe - Zj)]j,k=o,...,2n+N+i 

j>k>0 

holds for real zq, . . . , Z2 n +N+i- By taking the limit zq < z\ < . . . < z% n — > — oo and Z2n+N+i 
identity 



2n+N+l 

Yl S S n ( Z 3 - Z k) = Pf 
j>k>2n 



[Fjk]j,k=0,...,2n+N [9j]j=Q,...,2n+N 
[~ 9k\k=0,...,2n+N 



(2.50) 

-oo, we find another 
(2.51) 



where 



sgn(fc - j) (j, k = 0, . . . , 2n) 



Fjk - 



1 (j = 0, . . . , 2n; k = 2n + 1, . . . , In + N) 

-1 (j = 2n + 1, . . . ,2n + N; k = 0, . . . , 2n) 

sgn(z fc - Zj) (j, fe = 2n + 1, . . . , 2n + N) 

and & = g k = 1 (j, A; = 0, . . . , 2n + N). Substitution of Eq.(^5l]) into Eq.( |2.49| ) yields 

2n+N , 2n+N r[ „ , r , 

[-Fjklj,k=0,...,2n+N Wjl j=0,...,2ri+JV 

[ — Sfc]/c=0,...,2n+A' 



(2.52) 



p(zO,...,Z2n+iv) = J| yw(z 3 ) J| (Zj-Zfc)Pf 
j=0 j>k>0 



(2.53) 



The Pfaffian in the above can be represented as a quaternion determinant, due to the following theorem (which is 
essentially Theorem 2): 

Theorem 2' 

For even N, we can rewrite p(zq, . . . , Z2 U +n) as 



n+JV/2-l 

p(z , Z 2 u+n) = ( Y\. r i) s 2n+AT Tdet[/°™"(z i , 2fc)]j,fc=0,...,2n-HV- 



The quaternion elements are represented as 



S cv ™(z 3 ,z k ) P v ™(z 3 ,z k ) 
D e ™( Zj ,z k ) S cvca (z k , Zj ) 



(2.54) 



(2.55) 



and s, is defined in Eq.( p.22 ). The functions 5 cven , D cvcn and / cvcn are given in terms of 5, D and / defined in 
Eq.( [2.18 ), according to 



S CV011 (^,z fc ) = S(z,,z fc ) 
D even (z j ,z k )=D(z j ,z k ) 
I cvcn (z 3 ,z k )=I(z J ,z k ) 
Here # stands for a set of substitutions 



^2n+N{z k ) 



# S2n+N 
# 



(2.56) 



*2n+Af(Zj) - ®2n+N(z k ) 



# 



S2n+N 



R 3 {z) ^ R 3 (z) ^R 2n+N ( z ) (j = 0, . . . ,2n + N - 1). 

S2ri+W 



(2.57) 



It is straightforward to show that the quaternion element fj k cn (zj, z k ) in Theorem ^'satisfies the conditions imposed 
on fj k (zj, z k ) in Theorem 3. This means that we can inductively write 
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Rp(zo,..., z 2n + P ) 



r 3 )s 2n+N 



n 2 ™0 V W ( Z 3) Uflk>o( Z 3 - Z k) 



2n 



-Tdet[i-r n (^,^)] J ,fe=o,...,2n 



+p- 



Then the p- level correlation function ( 2.45 ) is written as 

S p (2 , ■ ■ ■ , Z 2n +p) _ T dct [/jfe Cn ( Z J' Z fe)]j,fc=0,...,2n+p 



p(x 1 ,...,x p ;{m}) 



We introduce a set of notations 



E (z , z 2n ) Tdet[/^ cn (zj, z k )]j,k=o,...,2n 



(2.58) 



(2.59) 



s 



IR 

jk 



S cvcn (z„z 2n+k ) (j = 0,. 



,2n) 



. , 2n; k = 1, 



S# = 5 cvcn (z 2n+J , * fc ) (j = 1, . . . , p; k = 0, . 

S RR = 5 cvc„ (z2r[+3; (j . fc = ^ jJ?) 

and similarly for D and /. A simplification 



-,p), 
,2n), 



(2.60) 



jk 



{-oo,z k ) = = S k , 

S 2n +N 

1 ^2n+N(z 2 n+k) 



S ' R = S cvca (-^,z 2n+k ) 



S 2n +N 



(2.61) 



T IR 

1 jk 



-1% = I c 



w x ®2n+N(z 2n +k) _ T R 

(.— 00, Z2n+/cj = = 4 , 

S 2n +N 



results from the definitions ( 2.18] ) and ( 2.56| ). Using Dyson's identity (A. 6) we can rewrite the correlation functions 



as 



Pf 



p(x 1 ,...,x p ;{m}) = (-l)^- i y 2 - 



Pf 



(-1) 



p(p-l)/2_ 







S n 


_j!R 




S IR ' 


(S H ) 


r 


D H 


-(S RI ) T 


D IR 


_jRI 




S RI 


_jRR 




S RR 


(S IR ) 


T 


D RI 


-(S RR ) 


T 


D RR 


Pf 




-I 11 


s 11 ' 










T D II 











s 1 


-I R 




S R 



-(5 7 ) T D" -{S RI ) T -{D RI ) T 
(I R ) T S RI ~I RR S RR 
-{S R ) T D RI ~(S RR ) T D RR 



Pf 



S 1 
-(S') T D 11 



(2.62) 



In the last line we have exploited a Pfafnan identity that holds for antisymmetric matrices A, B of odd ranks and a 
row vector v: 







V 






1 


Pf 


A 


V 


= Pf 


A 


1 




—v 1 ■ ■ ■ — V 1 


B 




-1 1 






Pf 



v 
-v T B 



(2.63) 



Likewise previous Subsection, we proceed to evaluate the component functions of the quaternion kernel fj k cn (zj, z k ) 
in the asymptitic limit where unfolded microscopic variables 



(2.64) 



V2Nz 


= Co-- 


= o, 








^2Nz 2j -i 


= C2i- 


-l = i/ij 


U = 


i,. 


.,n), 


V2Nz 2j 


= C2j 


= -Hh 


U = 


i,. 


. ,n), 


2Nz 2n+ j 


= 




(j = 


i,- 


■ ,P), 



11 



are kept fixed. We note that all elements of the sub-matrices that appear in the seco nd lin e of Eg . ( 2^62 ) are expressed 
as (derivatives or integrals of) analytic functions S(z, z') and &j(z) defined in Eqs.( 2.35 ) and ( 2.36 ), and ^j{z): 



1 d 



D ^ = 2dz- S{z ^ Zk) ' 



Sfk — S( X jl z k) + 



4> 



( X j) 



2a+N\ x j 
S2a+N 



S(-oo, z k ) - 



2a+N 



(Zk) 



S2a+N 



D 



BI 



Sjk — S(xj,Xk) + 



( x j) 



2a+N(Xj 
S2a+N 



S(-oc,x k ) 



2a+N 



(x k ) 



S2a+N 



S(— oo, x)dx — 1 



S(— oo, x)dx — 1 



1 d 



t rr _ o 
ijk - ~ z 



S(xj, x)dx — s,gw(xj — Xk) 



$2a+N(Xk) 



S2a+N 



S ( — oo, x)dx — 11 — 



4> 



2a+N 



( X j) 



S2a+N 



S ( — oo, x)dx — 1 



Utilizing an identity 



f°° i 
/ e^' 2 H 2k (u)du = 2 2k -^Y{k+ l -) 1 
Jo 1 



(2.65a) 
(2.65b) 
(2.65c) 
(2.65d) 
(2.65e) 

(2.65f) 
(2.66) 



an asymptotic formula 



fe-i 



-uV2 H Mdv _ 2 2fe - 1 r(fc) ^ , r(£ + i/2) 

H 2fc _i(«)d«- -j= ^(-1) r(£ + 1) 



< r(i + i/2) 
r(£ + i) 



2 2fc - 1 r(/c) f /¥_ (-i) fc 
2 ~ 2Vfc 



(*» i), 



and Eq. d2.38Q , we obtain Eqs.(|2.40|) (with 5, D, I in the second places replaced by S cvcn , D cvca , I cvcn ) and 

Si = S\ Z] ) ~ i — L= — cosCj, 



S"( Xj ) 



(-1) 



i+N/2 



■ COS A 



J ■ 



If = V2NI R ( Xj ) ~ (-l) n+Ar / 2+i y|sinA, 



(2.67) 

(2.68a) 
(2.68b) 
(2.68c) 



These matrix elements constitute the finite- volume partition function 



Z(M)^E ({-^=}) 



Pf 



const. 



S 1 
-(S') T D" 



(2.69) 



and the scaled spectral correlation functions 



p s (Ai 



p;{Al})= (72^)"^"-"A ;{ 7^ }) 



A,; 
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Pf 



(-1) 



p(p-l)/2. 



S 1 -I R S R 

-(s J ) T d" -{s RI ) T -cd ri ) t 
-(s R ) T r> RI -(s RR ) T r> RR 



Pf 



S J 
-(SO 7 D /7 



It satisfies a sequence 



p s ({A};/ii, . . .,//„, -/xi, . . . , -/u n ,0) ^i^ 00 p s ({A}; /ii, . . . ,fj, n -i, ••• j -Mn-1,0) 



/J. 71 _ i — - X. 



(2.70) 



(2.71) 



as each of the masses is decoupled by going to infinity. To illustrate this decoupling, we exhibit in FIG. 2 a plot of the 
spectral density p s (X; fi, — /i, 0) (p = l,n = 1). 

C. odd N f , odd N 

We finally consider the case with Nf = 2n + 1 flavors, {m} = (mi, . . . , m„, —mi, . . . , —m n , 0), and with odd TV. As 
mentioned in Introduction, this case is pathological because 



; { TO }) = -H p (xi, . . . , {m}) 



(2.72) 



and namely So({m}) = Z({m}) = 0. The quantities comput ed be lo w sho uld not be considered as correlation functions 
of a RME, but merely as multiple integrals defined by Eqs. ( 2. 44 )~(2. 46). 

As in previous Subsection, we define variables Zj by Eq.( 2.47 ), and denote the integrand in 2 p as 



2n+N 



. 2n+N 



2n+N 



p(z , z 2n +N) = H V w ( z ^ II ( z i _ Zfc ) IT s S n ( z i ~ Zfc )- 

j=0 j>k>0 j>k>2n 



(2.73) 



An identity 



2n+N 



YJ_ sgn(zj - z k ) = Pf[sgn(z fc - Zj)]j,fe=o,...,2n+iV 
i>fe>o 

holds for real zq, Z\, . . . , Z2n+./v- By taking the limit zo < z\ < . . . < z-i n — ► — oo, we find another identity 

2n+N 

Yl Sgn(Zj - Z k ) = Pf[-Fjfe]j,fc=o,...,2n+JV, 
i>fe>2n 



(2.74) 



(2.75) 



where F jk is defined in Eq. (|2.52D . Substitution of Eq. fl2.75D into Eq. fl2.73| ) yields 



2n+N 



, 2n+N 



p(z , . . . ,Z 2n +N) = Yl \jw{Zj) Yl (Zj ~ Z k )Pf[F j k]j,k=0,...,2n+N- 

j=0 j>k>0 



(2.76) 



The Pfaffian in the above can be represented as a quaternion determinant, due to the following theorem (which is 
essentially Theorem 1) 

Theorem 1 ' 



For odd N, we can rewrite p(zo, . . . , Z2u+n) as 



n+(JV+l)/2-l 

P(zq, . . . ,Z 2n+N ) = ( Yl r j)^ et [fjk(Zj,Z k )]j,k=0,...,2n+N- 

3=0 



(2.77) 
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The quaternion elements fjk{zj,Zk) are represented as 



fjk(Zj,Zk) 



S(Zj,z k ) I(Zj,z k ) 
D(zj,z k ) S{z k ,Zj) 



The functions S(zj,Zk), D(zj,Zk) and I(zj,Zk) are given by 

n+(N+l)/2-l 



S(zj,z k ) 
D(zj,z k ) 

I(Zj,Z k ) 



E 

1=0 
n+(N+l)/2-l 

E 

1=0 



^2t(z )^2l+l{zk) ~ ^2t+\{Zj)^2i{Zk) 



n 



n+(7V+l)/2-l 



E 

1=0 



$2l(Zj)$2l+l(*fc) - ^2l+\{Zj)^2t{Zk) 



(2.78) 



(2.79) 



As before, the above quaternion elements satisfy the conditions imposed on fjf.(zj,Zk) in Theorem 3, so that we 
can inductively write 



1p\Z0, ■ ■ ■ , Z2n+p) 



U p+(N+l)/2-l 
llj=0 



rii=o y/ w ( z i) Ujl k >a( z i ~ z k 



-Tdet[fj k {zj, Zk)]j,k=o,...,2- 



n+p- 



(2.80) 



In order to circumvent the vanishing of the partition function, we regard zq,zi,. . . , z 2n as generic variables, and 
define the 'p- level correlation function' 



p(xx, - ■ ■ ,x p ;{m}) 



Sp(z , • • • , z 2n+p ) _ Tdet[f jk (zj, z k )]j,k=o,- 



,2n+p 



E (zo, • • • , Z2n) Tdet[fjk(zj , Zk)]j,k=0,-,2r, 



(2.81) 



We introduce a set of notations 



Sjk = S ( Z ji Z k) 



(j,k = 0,..., 2n) 



Sjk = S(z j ,z 2n+k ) {j = 0, . . . , 2n; k = 1, 



,p), 



s jk = S(z 2n + 3 ,z k ) (j = 1, • • ■ ,p; k = 0, . . . , 2ra), 



Sjfe = S^n+j^n+fc) (j, fe = 1, • ■ ■ , p) 

and similarly for D and /. A simplification 

-lZ I = I IR (-OC,Z2n+k)=Ik H , 



jIR 



s 



kj 

IB. _ C IR 
jk 



(2.82) 



(2.83) 



S JK (-oc,z 2 „ +fe ) = Sf 



resu lts a gain from the definition ( 2.7S ). Using Dyson's identity ( A.6) we can rewrite the 'correlation functions' as 
Eq.(^62|). 

Likewise two previous Subsections, we proceed to evaluate the compo nent functions of the quaternion kernel 
fjk(zj,Zk) in the asymptitic limit where unfolded microscopic variables ( 2.64 ) are kept fix ed. W e note that all 
elements (other than S 1 , S R , and I R ) of the sub-matrices that appear i n the second line of Eq. (2.62) are expressed as 
(derivatives or integrals of) an analytic function S(z, z') defined in Eq.( 2.35| ) (with N replaced by N + 1), with real 
and/or imaginary arguments: 



D*l = \-^—S{zj,Zk), 



lk 2 dZi 



D fk = 2~dx~-^ <yXj,ZkS> '' 



(2.84a) 
(2.84b) 
(2.84c) 
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Sfk R — S(Xj,Xk), 



lf k R = 2 



S(xj, x)dx — sgn(xj — Xk)- 



We obtain Eqs.fl2.40p and 



(_l)«+(iV+l)/2 

S = S (zj) -== sin£j, 



S R ( Xj ) 



(_l)"+(A'+l)/2 



27T 



sin A ., 



If = V2NI R (x,) ~ (-l^+^+D/^cosAj 



(2.84d) 
(2.84e) 

(2.84f) 



(2.85a) 
(2.85b) 
(2.85c) 



Note the phase shifts of the trigonometric functions between Eq.( |2.85 ) and its even- TV counterpart (2.6S). These 
matrix elements constitute the 'finite-volume partition function' 



2({ M }) = So({- 



>2N 
Pf 



}) 



const 



S 1 
-(S') T D" 



and the 'scaled correlation functions' 

Pa 



5 (Ai,...,A p ;{/i}) = \Q=) P( 



(2.86) 



Pf 



(-1) 



p(p-l)/2. 



2N V2N V2N 

S 7 -I R 
-(S') T D /7 -(S flJ ) T 
(I fl ) T S RI — I RR 

-(s R ) T n RI -(s RR ) T 



s R 

-(T> RI ) T 

gRR 



Pf 



s 1 

-(S 7 ) T D" 



(2.87) 



We again remind the reader that under the identification (2.64), the above 'finite-volume partition function' vanishes 
and the 'correlation function' diverges as its Pfaffian denominator vanishes. 



III. SYMPLECTIC ENSEMBLE 

For (3 = 4, we concentrate on even flavor cases for a technical reason, and treat the following two cases separately: 

A : {to} = (mi, toi, . . . , m a , m a , —mi, —mi, . . . , —m a , —m a ), 

B : {to} = (mi, toi, . . . , m a , m a , —mi, —mi, . . . , —m a , —m a , 0, 0). 

The j3 = 4 case with an odd number of fermions will not be treated in this Article. 



A. N f = mod 4 

We first consider the case with Nf = 4a flavors and {to}=(toi,TOi,. . .,m a ,m a ,— mi,— mi,. . .,— m a ,— m a ). We 
express the partition function ( 1.12 ) of the RME in terms of eigenvalues {xj} of H: 



Z({m}) 



i fOO pOO " " / " 

h— 1H' ii e ~ 2x? n " + '"^ ) - ^) 4 

J -°° J -°°i=l 7=1 V fc=l 



(3.1) 



j>k 
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The p-level correlation function of the matrix H is defined as 

v 

p(xi,...,x p ;{m}) = (Y[tr8(xj - H)) 

3 = 1 

E p (x 1 ,...,x p ;{m}) 



~ ({m}) 



3 p (xi,...,x p ;{m}) = 



1 poo poo N N , a 

whvJ-J II n(- II 

PJ-J-OO J~°G j=p+1 k=1 



2 i 2 - 

Xj +m ki 



N 



Y[(xj ~x k y 
j>k 



(So = Z). We define new variables Zj as 

Z2j = —irrij 

Then the multiple integral ( |3 . 3| ) is expressed as 

^p(xi, ■ ■ ■ ,x p ; {z}) = — ^ - 2q - 
Uj=l W ( z 3)U.i>k( Z J ~ z w 



(j = l,...,a). 



(3.2) 
(3.3) 

(3.4) 



lj>k 

poo poo 2a+N 2a N 2a N 2a N 

-^t/ ■■■/ n ^•n^)n^) 2 iife-^)iiiife-^) 2 ii( 

PJ-J-oo J -co 0=2a+p+ i j-=l j= i j>k i=l Jfe=l j>fe 



(TV- 



where u>(z) = c z . 

Let us denote the integrand in Eq.([3.5[) as 



Xj X k J 



(3.5) 



2a N 2a N 2a N 

p(zi, . ■ .,z 2a ;xi, . . . ,x N ) = Y]_ w{zj) \\ w(xj) 2 \\_{zj - z k ) JJ J\(xj - z k ) 2 Y\_(xj - x k ) 4 

3 = 1 j=l j>k j=lk=l j>k 



(3.6) 



The above expressions can be represented as a quaternion determinant. In doing so, we introduce monic skew- 
orthogonal polynomials Qj(z) = z° + • • • that satisfy 



where 



{Q2j,Q2k+i)Q — -(Q2k+i,Q2j)g = qjSjk, 

(Q23, Q2k}Q = (Q23+I, Q2k+l)Q = 0, 



(f,9)Q= / dzw(z) 2 {f(z)gi(z)-f(z)g(z)). 



(3.7) 



(3.8) 



Explicit forms for the skew-orthogonal polynomials and their norms associated with the Gaussian weight w(z) are 
known in terms of the Hermite polynomials pOI]: 



2 3j + l/2 

Q^+l( Z ) = 2 3j+3/2 g 2j+l(V^z), 

q 3 = 2-^- 1 ' 2 ^{2j + 1)! . 
It can be readily seen that (cf. Eq.(26) of Ref. p9]) p(z\, . . . , z 2a \ xi, . ■ ■ ,Xn) is represented as a Pfafnan: 



(3.9) 



p(zi, . . .,Z2a',X!, . . .,x N ) = det 



" *oOi) 




• *2A r +2a 




^0(z2a) 


*l(«2a) ' 


' *2Af+2«- 




%(xi) 




*2JV+2a- 




*o(a:i) 


*i(a;i) • 


' ^2^+20- 


-loo 


%M 




• VI/' 




. ^o(x N ) 


*i(xjv) • 


' ^2N+2a- 


-i(ziv) . 
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a+N-1 

{ n ^) pf 

3=0 



-I(z 2 j,Z 2 k) I(z 2 j,Z 2 k-l) 
I{z 2 j-l 1 Z 2 k) -I(z 2j -i,z 2 k-i) 
-S(z 2k ,Zj) S(z 2k -i,Xj) 

I(Xj,Z 2k ) -I(Xj,Z3k-l) 



j,k=l, 



j=l,...,N;k=l, 



S(z 2j ,x k ) I(z 2j ,x k ) 
-S(z 2 j-l,X k ) -I(z 2j -i,x k ) 

D(xj,x k ) S(x k ,Xj) 
-S(xj,x k ) -I(xj,x k ) 



Here 



^ j (z) = w(z)Q j (z), 
and the functions S(x,y), D(x,y), I(x,y) are given by 



s( x y) = " v -1 

3=0 



1j 



D[x v) _ *^)*W*> : *W^M _ B_ S(X:A 

3=0 



dx 



I(x,y) = - 



a+N—l 

E 

3=0 



^ 2] {x)^ 2]+1 {y) - ^ 2j+1 (x)^ 2j (y) 



S(x, z)dz. 



Using Dyson's identity ( [A.q ), we can convert the Pfafhan into a quaternion determinant 

a+N-l 



p(zi,...,z 2a ;xi,...,XN) = ( { J ftjTdet 



3=0 



[g(z 2j -i, z 2j ; z 2k -u z 2k )] jtk=li ^ i0l [h(z2j-u Z2j-,Xk)} 
[h( z 2k-i,z2k;xj)] j=1 ^ N . k=li <a [f(xj,x k )} jik=lt _ 



where 



d{ z 2j-l, Z 2 j] Z 2k -l, Z 2 k) 
h(z 2 j-i,Z2j-,Xk) = 



-Ifoj-UZik) I{z 2 j-l,Z2k-l) 
-I(z 2 j,Z 2k ) I(z 2 j,Z 2k -l) 

S{z 2 j- 1 ,X k ) I(z2j-l,Xk) 

S(z 2 j,x k ) I(z 2j ,x k ) 

S{x ,x k ) I(xj,x k ) 
D(xj,x k ) S(x k ,Xj) 



f(Xj,X k ) 

and h stands for the dual of h (see Appendix A) . The skew-orthogonality relations ( |3.7[ ) lead to 

f(%\ x)f(x, x")dx = f(x\ x"), 
h(z, z'; x)f{x, x')dx — h(z, z'; x'), 
h(z, z , x)h(w, w'; x)dx = g(z, z'\ w, w'), 



which mean that the condition of Theorem 3 is satisfied. Therefore we can inductively write 



^p{x\ : • ■ • , Xp : Z\ . . . . , z 2a *j 

■w-ra+N-l 

n^o 13 

U 2 3 Zi w ( z j)U 2 j> k ( z 3 - z k 



-Tdet 



[g(z 2 j-l, Z 2 j] Z 2k -1, Z 2k )]j,k=l,...,a [h( z 2j-l, z 2 j',Xk)]j=l,...,a;k= 
[h(z 2 k-l, Z 2k ] Xj)]j=l,...,p;fc=l,...,a [f(Xj,Xk)]j,k=l,...,p 



Then the p-level correlation function (3.2) is written as 
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p{x 1 



^p{p^l ) * • • 5 *^P5 -2-1 9***3 ^2ri 



Tdct 



"0V^1> • • • i Z 2nj 

[g( z 2j-l, Z2j',Z2k-li %2k)]j,k=l,...,a [h[Z2j-l , ^2j 5 ^fe )jj=l,...,a;ft=l,...,p 

[/i.(22fe-i,^2fe;a;j)]j=i,...,p;fe=i,...,« [/fe'j arfc)]j,fc=i,...,p 



Tdet[3(z 2j -_i, z 2j ; z 2k -i, z 2 k)]j.k= 



We introduce a set of notations 



(3.17) 



(3.18) 



and similarly for D and /. Using Dyson's identity ( |A.6|) again, we can convert the quaternion determinant back to a 
Pfaffian, so that the correlation functions read 



Sj k = S{z j ,z k ) (j,k = l,...,2a), 

s jk = S(zj,x k ) (j = 1, ... ,2a; k = l,...,p), 

Sf k R = S{ Xj ,x k ) {j,k = l,...,p), 



Pf 



p(x 1 ,...,x p ;{m}) = (-ir^/ 2 - 



-I 11 -I IR S IR 
fjiR\T _jRR gRR 

-(S IR ) T -(S RR ) T D RR 



Pf[-J"] 



(3.19) 



Now we proceed to evaluate the component functions of the quaternion kernel, which are expressed as (derivatives 
or integrals of) the analytic function S(x,y) with real and/or imaginary arguments, in the asymptotic limit where 
unfolded microscopic variables 



V2Nz 2j -i = C2j-i = if-ij (j = l,...,a), 
V2Nz 2j = C2j = -i(J,j (j = 1, • • ■ , a), 
V2Nx 3 =\ 3 (j = l,...,p), 

are kept fixed. We note that S{x, y) has a compact expression 

e~ x2 ' y2 H 2a+2N (V2x)H 2a+2N - 1 (V2y) - H 2a+2N - 1 (\/2x)H 2a+2N {^y) 



(3.20) 



S(x,y) = 



2 2a + 2N y/^T(2a + 2N) 



H 2a + 2 N (V2y) 



x-y 
H 2a+2N -i(V2u)di 



(3.21) 



In the second line of Eq.(3.21), we have singled out the unitary scalar ker nel an d applied to it the Christoffel-Darboux 
formula. Substituting asymptotic formulas for the Hermite polynomials ( [2.38 ), one can show that S(x,y) approaches 
the sine kernel fief : 



1 ,S(- C ^ 



'2N 



<2N V2N 



2n((-C) -^( 2 (C-C))- 



(3.22) 



The Pfaffian elements in the the asymptotic limit, after taking into account an unfolding by the factor \ 2N, are then 
expressed in terms of K(C,): 



T 77 


= I( Zj ,z k )~ /V(2C)dC, 
Jo 




(3.23a) 


alR 




-A*)), 


(3.23b) 


rlR 
L jk 


= I(z j>Xk )~ [ V(2C)dC, 
Jo 




(3.23c) 


a RR 
=>jk 


^^S(x J ,x k )^K(2(X J 


-A*)), 


(3.23d) 




= ^D(x j ,x k )~2K'(2(\ j 


-A fc )), 


(3.23e) 


-iRR 
L jk 


= I(xj,x k )~ [ J K{2X)dX. 
Jo 




(3.23f) 



18 



These matrix elements constitute the finite- volume partition function 



Z({m}) = E ({- 



'2N 



}) 



= const. 



Pf[-I u ] 



and the scaled spectral correlation functions 



p s (Ai 



A„iM)=(ii=)'p(-F= -^=;{ 



2N 
Pf 



(-1) 



p(p-l)/2_ 



/2JV 
-I" 



:}) 



/2iV 

gRR 

-(s IR ) T -(s RR ) T n RR 



Pf [-!"] 



(3.24) 



(3.25) 



In the quenched limit /ii, . 



oo when the ratio of two Pfaffians is replaced by a minor Pf 



—i s 

_( S HH)T DBS 



the correlation functions approach those of the Gaussian symplectic ensemble [}40| 
sequence 

,o s ({A}; • • • , Mn, Mn: ~Mi j _ Mii • • • j ~Mn) 

Ps({A}; /ill Mli ■ ■ ■ 7 Mn-l, Mn-l) — Mil — Ml> ■ ■ ■ 7 — Mn-1) — Mn-l) 



By the same token, it satisfies a 

(3.26) 



as each of the masses is decoupled by going to infinity. To illustrate this decoupling, we exhibit in FIG. 3 a plot of the 
spectral density /0 S (A; /i, /i, — — //) (p = 1, a = 1). 



B. N f = 2 mod 4 



Next we consider the case with Nf = 4 a + 2 flavors and {to} = (mi, mi, . . . , m Q , m a , mi, —mi, . . ., — m a , — m a , 0, 0). 
We express the partition function ( 1.12j ) of the RME in terms of eigenvalues {xj} of H 



N 



N 



The p-level correlation function of the matrix _ff is defined as 

p 



P ;{m}) = (l[tTd(x j -H)) 

3=1 

_ a p (xi,...,a:p;{m}) 



S ({m}) 



H p (a;i,...,Xp;{m}) = 



OO /"DO 



(N-p)\ 

(Sq = Z). We define new variables Zj as 



oo J — oo 



AT JV 

n 

' 3=P+1 J=l 



fe=l 



iY 



j>fc 



2 = 0, 

■22.7-1 = imj 
22 j = —irrij 



(j = l,...,a). 



Then the multiple integral ( 3.29 ) is expressed as 

^pi^ll • ■ • : Xp : {z}) 



n 2 = n, 2 >fc> ( z J - **) 



(3.27) 



(3.28) 
(3.29) 



(3.30) 



,.00 <.oo 2a+W 2a JV 2a N 2a N 

v — t / • • • / n ^n^n^) 2 n (^-^nn^-^) 2 ]!^-^) 4 ' ( 3 - 31 

l Pj- J-oo J-oo j=2 q+P+1 .7=0 j"=l j>fe>0 j=l fc=0 j>k 
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where w(z) — e z . 

Let us denote the integrand in Eq. ( |3.31 ) as 



N 



N 2a 



N 



p(z , . . .,z 2a ;x 1: ■ ■ -,x N ) = Y[w(zj) Y[w(xj) 2 Y[ ( z j - z k) J! Y[( x i ~ z k) 2 

j=0 j = l j>k>Q j=l k=0 j>k 

It can be readily seen that p(zq, . . . , Z2 a ', x±, . . . , xn) is represented as a Pfafhan: 



p(z , . . .,Z2a,Xi, . . .,x N ) = det 



^n(zi) *i(zi) 

^o(z 2a ) Vi{z 2a ) 
*o(zi) 



%(x N ) %{x N ) 



^2N+2a(zo) 



^2N+2a(z2a) 

%N+2 a ( X l) 
^2N+2a(xl) 



a+N-1 

( II <U)W 

3=0 



[Ak]k=l,...,a 
[Hkj]j=l,...,N;k=l,...,a 



%N+2a( x N) 
^2N+2a(xN) 

[Aj]j=l,...,a [Hjk]j = l,...,a;k=l,...,N 

£1 [Bk]k=l,...,N 
[Bj]j=l,...,N [Fjk]j,k=l,...,N 



where 



G jk 
F jk 
Bk 

n 



-I(z 2 j,Z 2 k) I{z2j,Z2k-l) 
I(z 2 j-l,Z2k) -I{Z2j-l,Z2k-l) 

S{z 2] ,x k ) I(z 2 j,x k ) 

-S(z 2 j-l,X k ) -I(Z2j-l,Xk) 

D(xj,x k ) S{x kl Xj) 
-S( X j,Xk) -I(xj,x k ) 

^2N+2a(z 2 j) I(z 2j ,Z ) 
-^2N+2ct{Z2j-\) -I{Z2j-l,Zo) 

-^2N+2a( X k) ^2N+2a(x k ) 

-S(z ,Xk) -I(z ,x k ) 

^ 2 N+2a(zo) 
-^2N+2a(z ) 



The functions S, D, I are given in terms of 5, D and / defined in Eq.(3.12) according to 

^2a+2N{x) 



S(x,y)=S(x,y) 
D(x,y) = D(x,y) 
I(x,y)=I(x,y) 



S2a+2N 



S2a+2N 



and 



Here * stands for a set of substitutions 



h( z ) ^ Qj( z ) 



j 



S2a+2N 



9j(x)dx. 



,+2n(z) (j = 0,...,2a + 2N-l). 
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Using Dyson's identity ( |A . q ) , we find 
p(z , . . . ,z 2a ;xi, . . . ,x N ) 

[9( z 2j-l, Z2j', Z 2 k-1, Z2k)]j,k=l,...,a [d(z2j-l , Z2j)]j=l,...,a [h( z 2j-l, Z2j ', %k )] j = l,. . . ,a;k 

[a(Z2k-l,Z2k)]k=l,...,a & [b(x k )]k=l,...,N 

[h( z 2k-l, z 2k', x j)]j = l,...,N;k=l,...,u [H X j)\j=l, — ,N [f( x j i x k)]j,k=l,...,N 



a+N-1 

( n ^ m " 1 

3=0 



where 



g(z 2 j-l, Z2j] Z 2 k-1, Z 2 k) 

h(z2j-i, z 2 j]Xk) 
f(xj,x k ) 

a(Z2j-l,Z2j) 

b(xk) 



-I(z2j-l,Z2k) I{z2j-l,Z2k-l) 
-I(z 2 j,Z2k) I{z2 3 ,Z 2 k-l) 

S(Z2j-l,Xk) I(Z2j-l,Xk) 
S(z 2 j 1 X k ) I(Z2j,X k ) 

S(xj,Xk) I{xj,x k ) 
D(xj,x k ) S(x k ,Xj) 

^2N+2a(z2j-l) I(z2 3 -l,Z a ) 

S(z ,x k ) I(z ,x k ) 

-%N+2a( X k) ^2N+2a{ x k) 
^2N+2a(z ) 

^2N+2a(z ) 



and a, b, h stand for the duals of a, b, h. The skew-orthogonality relations ( |3.7| ) lead to 

f(x',x)f(x,x")dx = f(x',x"), 

b(x)f{x,x')dx — b(x'), 
h(z, z'; x)f(x, x')dx — h(z, z'; x') 1 
b{x)b(x)dx = u), 
h(z, z'; x)b(x)dx — a{z, z'), 
h(z, z', x)h(w, w'; x)dx — g(z, z'; w, w'), 



na+N-l 
3=0 1j 



which mean that the condition of Theorem 3 is satisfied. Therefore we can inductively write 

^■p \XX j • • • ) Xp 5 Zq , . . . , Z2a ) 



xTdet 



]lj=0 W (Zj) Uj>k> ( Z 3 - Z k) 

[g( z 2j-l, Z2j', Z2k-1, Z2k)]j,k=l,...,a [ a { z 2j-l , Z2j)]j=l,...,a [h( z 2j-l, Z2j', Xk)]j=l,...,a;k=l,...,p 

[a(z2k-l,Z2k)]k=l,...,a U [b(xk))k=l,...,p 

[h(z 2 k-l, Z2k] 3Cj)]i=l,...,p;fc=l,...,a ${ x j)]j=l,—,P [f(xj,X k )}j,k=l,...,p 



Then the p-level correlation function (3.28) is written as 

^p{X\ , . . . , Xp] Zq, . . . , Z2n) 



p(xi,...,x p ;{m}) 



So(zo) ■ • • 7 z 2n ) 



Tdet 


[g( z 2j-l, Z2j', Z2k-1, Z2k)]j,k=l,...,a [a.(z 2 j~l , Z 2 j )] j=l,. . . >a [/l(z 2 j_l, Z 2 j] X k ) 
[a(Z2k-l,Z2k)]k=l,...,a u [b(x k )]k 
[h(z2k-l,Z2k;%j)]j=l,...,p;k=l,...,a [K x j)]j=h-,P [f( x j,Xk) 


]j — l,...,a;k— l,...,p 
— l....,p 

j,k=l,...,p 


Tdet 


[g(z 2 j-l, Z2j] Z 2 k-1, Z2k)]j,k=l,...,a [a{z2j-l, Z2j)] j=l,...,a 
[a(z 2 k-l,Z2k)]k=l,...,a W 
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Now we make replacement of the elements of the quaternion kernel back to those defined in Eq.(3.12): 

S(x, y) -> S(x, y), D(x, y) -> D(x, y), I(x, y) -> I(x, y), (3.43) 

which does not change the values of the quaternion determinants in Eqs.( 3.41 ) and ( 3.42| ). We introduce a set of 
notations 



s jk = S(zj,z k ) (j, k = 0, . . . ,2a), 
S IR = S{ Zj ,x k ) (i = 0,...,2a; k = 1, 
Sil R = S{xj,x k ) (j,k = l,...,p), 



,p), 



(3.44) 



and similarly for D and /, and 



Ff 





U 


= o,. 


.,2a), 


S2N+2a 






U 


= 1,. 


• ,P), 


$2N+2u 




^2N+2a( X j) 


U 


= 1,. 




S2N+2a 









(3.45) 
(3.46) 



Using Dyson's identity (A. 6), we can convert the quaternion determinant back to a Pfaffian, so that the correlation 
functions read 



Pf 



p(x 1 ,...,x p] {m}) = (-iy^ 2 - 



-I 11 Q 1 —i IR s IR 

-(Q T ) T o ~{Q R ) T ~(p r ) t 

(I IR ) T Q R -I RR S RR 

-{S IR ) T P R ~{S RR ) T D RR 



Pf 



-I 11 Q 1 

-(QT o 



(3.47) 



Now we proceed to evaluate the component functions in the asymptitic limit where unfolded microscopic variables 



(3.48) 



V2Nz 


= Co = o, 










r 2Nz 2j -i 


= Oy-i = 




U = 


i,- 


. ,a), 


V2Nz 2j 


= &j = 




(? = 




. ,a), 


V2Nxj 






(j = 


i,. 


■ ,P), 



are kept fixed. We obtain Eq.(p.23|) and 



T}R 



1 



'2N 



'2N 



P 



R 



0. 



These matrix elements constitute the finite-volume partition function 



(3.49a) 

(3.49b) 
(3.49c) 



Z({ M }) = So({- 



/2AT 
Pf 



}) 



const. 



-I" Q 7 

-(QT o 



2\2 ' 



(3.50) 



and the scaled spectral correlation functions 
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p s (Ai, . . . , A p ; {n}) = \q= ) P( 



Ai 

/2N 



A,, 



{- 



M 



Pf 



(-1) 



p(p-l)/2_ 



V2iV V2iV 
-I 77 Q 7 -I 77? 
-(Q') T o -(Q*) T 
(I 77i ) T Q 7 * 
-(S IR ) T P R -(S RR ) T 



S IR 

-(p fl ) T 



Pf 



I 77 Q 7 

-(QT o 



(3.51) 



It satisfies a sequence 



/9 a ({A};/ii,/xi, . 

p s ({A};Mi)Mi) • 



• , Mm M") — Ml) — Ml) ■ • • J 
! Mn-1) Mra-1) — Ml) — Ml) 



-Mn,-Mn) 0,0) ► 

• .,-Mn-l>-Mn-l)0,0) 



(3.52) 



as each of the masses is decoupled by going to infinity. To illustrate this decoupling, we exhibit in FIG. 4 a plot of the 
spectral density p s (A; m, M) — M, — M> 0j 0) (p = 1, a = 1). 
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APPENDIX A: QUATERNION DETERMINANT 

A quaternion is defined as a linear combination of four 
basic units {1, e±, e%, e^\: 



q = <7o + q - e = q Q + q x e\ + q 2 e 2 + q^. 



(A.l) 



Here the coefficients qo , qi , q 2 , and q^ are real or complex 
numbers. The first part qo is called the scalar part of 
q. The quaternion basic units satisfy the multiplication 
laws 



1-1 = 1, 1 ■ e,- 



l = e s (j = l,2,3), 



e 2 = e 3 = eie 2 e3 = — 1. 



(A.2) 



The multiplication is associative and in general not com- 
mutative. The dual q of a quaternion q is defined as 



q = qo - q • e. 



(A.3) 



For a selfdual N x N matrix Q with quaternion elements 
qjk has a dual matrix Q = [cjkj]- The quaternion units 
can be represented as 2 x 2 matrices 



(A.4) 



We define a quaternion determinant Tdet of a selfdual Q 
(i.e., Q — Q) as 



"10" 




'0 -l" 


1 


, ei -> 


1 





" -i~ 




' i " 


62 -> 


-i 


, e 3 -> 


-i 



TdetQ= J] (-1)^II(«„6 



•g dQ )o, (A.5) 



Pes* 



where P denotes any permutation of the indices 
{!,..., N) consisting of I exclusive cycles of the form 
{a — >6— >c— > d — * a) and (— l) N ~ l is the parity 
of P. The subscript means that the scalar part of the 
product is taken over each cycle. Note that a quaternion 
determinant of a selfdual quaternion matrix is always a 
scalar. The quaternion determinant of Q can as well be 
represented by its 2N x 2N complex matrix representa- 
tion C{Q) §§: 



TdetQ = Pf [JC(Q)], J = l 



N 



1 

-1 



(A.6) 
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We adopt Cartan's convention on Riemannian symmet- 
ric spaces |l^| and supplement subscripts to indicate the 
ranks of Lie group constituents. 

We tacitly assume that a regularization employed in con- 



voluting the operators e and K in Eq.(1.21 



) does not in- 
volve non-universality. For a detailed discussion, see Ref. 
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FIG. 1. The scaled spectral density p s (C; Mi — m) for the orthogonal ensemble with two flavors (n = 1). 
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